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We investigated theoretically the Coulomb drag effect in coupled 2D electron gases in a wide 
interval of magnetic field and temperature 1/r Cw c C Bf/S, T < Bp, t being intralayer scattering 
time, ui c being the cyclotron frequency. We show that the quantization of the electron spectrum 
leads to rich parametric dependences of drag transresistance on temperature and magnetic field. 
This is in contrast to usual resistance. New small energy scales are found to cut typical excitation 
energies to values lower than temperature. This may lead to a linear temperature dependence of 
transresistance even in a relatively weak magnetic field and can explain some recent experimental 
data. 

We present a novel mechanism of Coulomb drag when the current in the active layer causes a 
magnetoplasmon wind and the magnetoplasmons are absorbed by the electrons of the passive layer 
providing a momentum transfer. We derived general relations that describe the drag as a result of 
resonant tunneling of magnetoplasmons. 

PACS numbers:73.20.Mf, 73.61.-r, 73.23.Ps 



I. INTRODUCTION 

When two 2-dimensional electron systems are placed in close proximity, then even in the absence of electron 
tunneling between layers the current in one layer (the active layer) will cause the current in the other (the passive 
layer). El This phenomenon is known as a frictional drag and is due to the interlayer Coulomb interaction which causes 
a momentum transfer from one layer to the other. If no current is allowed in the passive layer, the potential difference 
developes there to compensate the frictional interlayer force. The transresistance is measured as the ratio between 
the electric field developed in the passive layer to the current density in the active layer. This drag is of fundamental 
interest because it can be used as a sensitive probe of the screened interlayer Coulomb interaction and the form of 
the irreducible polarizability function x(w,q) within the layer. n 

Until recently, the experiments on frictional drag have been done in a zero magnetic field.EJ In this regime, the 
theoretical description is well elaboratedBH. The T 2 dependence of transresistance was explained by the phase-space 
arguments. mm 

A number of experiments on Coulomb drag in magnetic field have been also reported™. Most of them are done 
in Quantum Hall regime (T <C huj c ~ Ep , T, uj c being temperature and cyclotron frequency respectively) when the 
screening of the interaction and the polarizability function are determined by the states of the lowermost Landau 
level. The experiments have shown a strong filling factor dependence of drag resistance. It was strongly enhanced 
(compared to the H = case) when the Fermi level lay within the Landau level and was strongly supressed when the 
Fermi level was between the Landau levels. The recent theoretical workQ'Q (which treats the case T < A <C hu> c ~ Ep 
, where A is the width of the Landau level due to disorder) predicts a twin-peak structure of the transresistance as 
5.^ ■ a function of the magnetic field. This is due to the interplay between the screened interlayer interaction and the 
& ' phase-space available for the interlayer e-e scattering. 

In the present work, we investigate the drag effect in intermediate magnetic fields 1/t <C lu c *C Ep /h. In contrast 
to Quantum Hall regime, the e lectrons occupy many Landau levels. Still these Landau levels remain well resolved 
since their width A ~ hy/uj c /T <C hto c is much smaller than the level spacing. 

The electron density of states is strongly distorted in comparison with its zero-field value. That is why even in the 
"classical" regime Ep >T> hu) c and R c ^> d (R c = vf/uj c being the cyclotron radius, d being the distance between 
the layers) the polarization function Imx(a; ~ T, q ~ 1/d) which is responsible for the absorbtion of energy, differs 
strongly from its zero-field form. As a function of frequency it consists of a series of well resolved peaks at multiples 
of cyclotron frequency. We will see that this circumstance leads to rich parametric dependences of transresistance 
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on temperature and magnetic field. This is to be contrasted with the usual intralayer resistance that exhibits no 
anomalies except strongly supressed Shubnikov-de Haas oscillations and does not manifest the electron density of 
states. 

In the intermediate magnetic field at temperatures T > ?ko c weakly damped boson excitations become important. 
Those are magnetoplasmons with energies close to multiples of cyclotron energy. These excitations provide a new 
mechanism of Coulomb drag in the system. The momentum transfer is provided by magnetoplasmons being excited 
in one layer and absorbed into the other. We have found general relations that allow to present the magnetoplasmon 
contribution to the drag resistance as the result of resonant tunneling of magnetoplasmons. n 

The problem we consider in this paper is interesting also because the experimental dependences of transresistance 
in the indicated parameter region has a universal form cx B 2 T . This temperature dependence is of a highly nontrivial 
form because the common arguments based on consideration of the phase-space available for the scattering give the 
T 2 dependence. We show that new small energy scales (for example, A < T) play the role of characteristic cut-off 
energies which leads to a linear temperature dependence of transresistance in a relatively weak magnetic field. 

A bird's eye view of our results is provided in figure 1, where we present temperature and magnetic field exponents 
of the transresistance in five distinct parameter regions. We obtain detailed analythical results for the first three 
regions that correspond to T> hu c . In addition to slow T— and B— dependence the transresistance in the regions 
IV- VII (T <§; huj c ) exhibits oscillatory dependence on inverse magnetic field (Shubnikov-de Haas oscillations). The 
oscillations can be hardly investigated analythically except simplest cases. In the present paper we only present 
analythical results for the regions VI and IV and give estimations of the transresistance in the regions V and VII for 
typical filling factors. 

The outline of the paper is as follows. Irusfiction 2 we list the theoretical assumptions we made and present the 
method which is essentially the same as in&ffiffl. Details of polarization function and magnetoplasmon spectrum are 
presented in section 3. In section 4 we give a detailed description of magnetoplasmon mechanism of the drag. A 
phenomenological description of the resonant tunneling of magnetoplasmons is elaborated in section 5. We list our 
analythical results for high temperatures in section 6. Section 7 is devoted to evaluation of the transresistance at low 
temperatures. 



II. METHOD 



In the present paper, we cover parameter region A <C hto c <C Eq = Tivp/d. Here A stands for the width of the 
Landau level. It determines the maximun one-particle density of states in a magnetic field. The ratio between fiui c 
and T can be arbitrary. n 

We assume here that Landau levels acquire width due to scattering by impurities and, following^, treat the effect 
in the selfconsistent Born approximation (SCBA). This approximation is known to lift the difficulties related to the 
high Landau level degeneracy. In this approach, A 2 = (2/ir)fruj c h/T. This expression for A is valid for a short range 
((5-correlated) random potential. p. 

We also assume Coulomb mechanism of the draga, so that the d.c. drag current results from the rectification 
by the passive layer of the a.c. fluctuating electric field created by the active one. In diagrammatic language, 
the transconductance is given by a diagramm composed of thre&jjbody correlation functions, those are connected by 
Coulomb interaction lines (photon propagators) . It was argued incHj that under very general conditions the three-body 
correlation functions can be expressed in terms of electron polarization functions x(w, q)i,2 in each layer. n 

This yeilds the following expression for the diagonal element of the transresistivity tensor( Eq. 28 of Refill): 

xx _ fi 2 1 f d 2 q 2 f°° dw_ Vi 2 (q) |2 ImXi (g, gjlmXgQj, g) m 
Pl2 ~ 2e 2 mn 2 T J {2^Y q J 2tt 1 £(cu,q) 1 sinh 2 (^/2T) ' [> 

Here ni, n 2 are electron concentrations in the layers, Vi 2 (q) is the Fourier component of interlayer Coulomb interaction 
and £(w,q) describes (dynamical) screening of this interaction. Electrostatics gives V\ 2 = V(q) exp(— qd), V(q) = 
2ire 2 /eq ,e being the bulk dielectric constant, and 

E(u, q) = (l + V(q) X l)(l + V(q) X2 ) - V 2 2 {q)XiX2- (2) 

It has been arguedB that Eq.l is valid for arbitrary magnetic field provided q^ 1 -C /,i? c , I being mean free path. 
Our checks confirm that, so we use the Eq. |l| in our calculations in intermediate magnetic field regime. 

As a inference, we give here the expression for transresistance in the absence of magnetic field for temperature 
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_ C(3) h 1 T 2 a 2 B r . 

P12 — o J4 72 ' 

Here £(3) ~ 1.202 is the Riemann zeta- function, Vfx,Vf2 are the corresponding Fermi velocities and as is the Bohr 
radius. 

To make use of Eq. [l], we shall evaluate the polarization function x(w, q; B). That we do in the following section. 

We will assume that the layers are macroscopically identical. We keep the indices 1, 2 that label the layers in the 
formulas solely for the sake of physical clarity. The exception will be the discussion of Shubnikov- de Haas oscillations 
in the regions IV and VI. We allow there for different filling factors in the layers. 

III. POLARIZATION FUNCTION AND MAGNETOPLASMONS 

We start with the following expression for the imaginary part of the polarization function: 

/+°o ^ £ 
— [n F {e) -n F {e + u))]lmG r n (e)lmG r m {e + w), (4) 
-°° 7r 

Here v — mo/irh 2 is the 2D thermodynamic density of states in the absence of magnetic field, np is the Fermi 
distribution function, and G r n is the retarded Green's function of the electrons in the nth Landau level. In the above 
formula, we have taken into account that under the conditions considered in this paper (T,Klo c -C Ep) only big 
Landau level numbers are important. Thus the bare vertex function is reduced to its quasiclassical form, which is 
Bessel function of the argument qR c . In the limit A — > (no disorder) this expression is equivalent to semiclassical 
approximations employed iroEij. 

The expression (0) disregards vertex corrections due to disorder. This is safe since we always assume that VFq 3> 1/r 
and qR c 3> 1. We will also disregard the rapidly oscillating part of Bessel function squares at qR c 3> 1, that is, we 
assume that J^iqRc) ~ ^/^qR c - We discuss the relevance of this assumption in Appendix. 

Using the SCBA expression for ImGQ in the limit of big n, 



^( e ) = -|^/i-(^e(i-(^) 



(•5) 



where is the step-function and e„ = (n + l/2)?kj c ,we obtain from Eq.(|J) that 



knx(w,q) = v-^rr — Y] [nF(en) - n F {e„ + w)] X im {— — ^" + huj y ( 6 ) 
7r z A qvp * — ' 2 A 



Here we define dimensionless function Xi m (x) = (4/3)[(l + x 2 )£ , (Vl — x 2 ) — 2x 2 F(\/l — x 2 )] for \x\ < 1 andX !m (a;) — 
otherwise. The functions F(x),E(x) are the complete elliptic integrals of the first and second kind, respectively. 
Note that J_ x dxX lm (x) = tt 2 /8. The expression for Im\(uJ, q) assumes different forms depending on the temperature. 
In the most interesting case T 3> huj c and qvF 3> uj (the latter inequality is equivalent to T -C Eq, because the 
characteristic q ~ 1/d and the frequency cannot be larger than the temperature) we get from Eq. (||) 

Imx(w,g) = v-A-— V X im Q. (7) 
7r z A qvp * — ' 

j=-co 

Since xi 10 ) is an analytical function of uj, we can easily obtain the real part from Eq.(Q). In the vicinity of the j-th 
cyclotron resonance, \u> — ju> c \ <C uj c , it reads: 

t> / \ . 4w c w u~juj c 1 f +1 dyX im (y) 

Rex{^,q) = v + v^— X r& (— - — ); X re {x) = -v.p. / — — . 8 

7r^A qvp 2A 7r 7-i [V ~ x ) 

The functions X re (x) and Xi m (x) are plotted in Fig. 2. 

In the opposite limit of A -c T -C hu> c and u < 2 A; qR c > 1 we obtain from Eq.(||): 

lmx(ui,q) = "^f^^i^-^-M 1 - /«)' ( 9 ) 
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where /„ = 1/(1 + exp(e„ — A*)/T) being the filling factor of the n-th Landau level in the layer, being chemical 
potential. 

Finally, at T < A we can set n = m and integrate expression (^) over e in close vicinity of [i. This gives 

Im *=^X2 1 ^2 ( 10 ) 

ir z A z qvp A z 

Although the magnetoplasmon modes of a 2DEG have been extensively studied0, very little attention has been 
paid to their properties at high frequencies in short-wave limit. Since those are of interest for us, we investigated 
them in some detail. The dispersion curves of the magnetoplasmon modes in the case of the weak damping are 
determined by the equation Re£ = 0, £ being given by Eq.||. From this equation we get {Re\/v) + z± = 0, 
where z± = qaB/2(l ± exp(— qd)) <C 1. Using Eq.(||) for Re\, for each j we obtain two solutions corresponding to 
two (±) magnetoplasmon modes. First thing to note is that under our assumptions qas <C 1 so that on a large 
frequency scale za-, can be safely omitted. The resulting equation Rex = determines a series of Burstein-like 
magnetoplasmonsEJ with frequencies that in the interesting region of q are close to cyclotron harmonics juj c - (see Fig. 
3) If q < qo = — 4j X re (l)uj 2 /ir 2 Avf — 0-I9joj c / R C A the root of dispersion relation lies beyond the electron adsorption 
bands and the magnetoplasmon is not damped. These magnetoplasmons are of no interest for us since they do not 
talk to electrons and thus cannot participate in drag. If q 3> qo magnetoplasmons lie deep in the electron adsorption 
band and are strongly damped, so that their conribution to drag cannot be distinguished from the contribution of 
electron-electron scattering. This is why we concentrate now on a close vicinity of go (right panel of Fig. 3). 

In this vicinity x( w j q) can be expanded in Taylor series in terms of w = lo — jw c — 2A, n = q — qo, assuming that 
w <C A, n <C qo, 

x/^c^ + ^ + ic^yei-,) (ii) 

Here are numerical constants characterizing behaviour of Xi m ,X re near x = 1, C% ~ 6.23, C2 — 2.19. This 
determines the dispersion law of magnetoplasmons 

2A 

w± = —?r(z± + K/qo) (12) 
(^2 



and their damping 



r (13) 



where z is taken at q = qo- Symmetric and asymmetric modes are split by 



5. = A-^—. (14) 
02 smhqo« 



We see that r <§; w <C A, Sui <C A, Sui can be comparable with T. 



IV. MAGNETOPLASMON CONTRIBUTION. 



In this section we consider the magnetoplasmon mechanism of Coulomb drag. 

In the absence of magnetic field there are two plasmon modes in double-layer system, the one with the electron 
densities in the two layers oscillating in phase (the optic mode), and the other one where the oscillations are out 
of phase (the acoustic mode)E3. It was pointed out inB that the drag effect can be greatly enhanced by dynamical 
" antiscrecning" of the interlayer interaction due to coupled plasmon modes. Since the plasmon modes lie beyond 
the T = particle-hole continuum, temperatures of the order of Fermi energy are required for a large plasmon 
enhancement of the drag effect. Only then the thermally excited electrons and holes with plasmon velocities provide 
sufficient damping of the plasmon modes and thus facilitate plasmon interaction with electrons! 

In the revised situation of intermediate magnetic field, the magnetoplasmons have even better chances to enhance 
the drag. First, there are many modes and their typical energies are of the order of Hlu c . Therefore these modes can be 
excited at temperatures much lower than Fermi energy. Second, the magnetoplasmons in our model acquire natural 
damping: due to finite Landau level width they may lay within the particle-hole continium. The finite temperature 
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without disorder does not lead to magnetoplasmon damping and to the drag effect. This is in contrast to the situation 
without magnetic field, where lm\ at the plasmon frequency was calculated for collisionless plasmaa. 

The magnetoplasmon mechanism of the Coulomb drag, when the current in the active layer causes a magneto- 
plasmon wind and the magnetoplasmons are absorbed by the electrons of the passive layer leading to transfer of the 
momentum, must be quite general one. In this section, we evaluate magnetoplasmon contribution with using Eq. ^. 
It turns out that the answer can be expressed through only two quantities for each double plasmon mode: frequency 
splitting Suj and damping T. Any concrete model would only set specific expressions for T and 5u>. This clarifies the 
physical meaning of Eq. |lj. 

To prove this, we rederive the result in the next section in a phenomenological framework. 

Let us expand xt^i l) around the frequency uj(q) at which Rex = 0, 

X(w,q) =x'w + ix", (15) 

w being frequency deviation. Here x' — R e dx/duj {&{<})), x" = I m x( w ( ( 7))- Expression for £ reduces to the form 

£ = (V 2 (q) - V? 2 {q))(wx' + ix" + vz+Xwyd + i X " + vz-) (16) 

Consequently, the integrand in Eq. |l| has a sharp maximum near ui{q) as a function of ui. This suggests we can now 
integrate over w in infinite limits. 

Eq. [l^ determines magnetoplasmon spectrum and suggests that mode splitting Suj = v{z- — z + )/x', damping 
r = 2x"/x'- All factors V12, V can be absorbed into these two quantities. Typical w that contribute to the integral 
are of the order max(r, Suj). The approximation is valid if x and statistical factor sinh 2 (u;/2T) do not change much 
in this frequency window, that is ui(q) ^> Su, dT/duj <C 1, max(r, Suj) <C T/h. 

Provided these conditions are fulfilled, we can reduce the expression for the transresistance to the following elegant 
form: 



XX 

P12 



d 2 qq 2 V 1 (M 2 r 



16e 2 nm 2 T J (27r) 2 J-f sinh 2 (^(g)/2T) [(fc) 2 + T 2 ] ' 



Summation over modes means summation over all possible roots of Rex = 0. 



V. RESONANT TUNNELING OF MAGNETOPLASMONS 



We give here another derivation of this formula which clarifies its physical meaning. To describe resonant tunneling 
of plasmons between the layers, we introduce for each plasmon mode a density matrix pij =< b\bj >. Here i,j = 1,2 
label the layers, b\ b are boson creation/annihilation operators. In the absence of dissipation, that is, plasmon emission 
and absorption, the density matrix obeys the equation 

^=iJ2(H il Pi j -p U H lj ) (18) 

For identical layers, diagonal elements of the Hamiltonian are equal to each other and dissappear from the equation. 
The non-diagonal element that is responsible for plasmon tunneling between layers can be readily express in terms of 
splitting 5u> between symmetric and assymmetric plasmon state: H12 = H21 — Slo/2. 

The dissipation takes place independently in each of the layers. It contributes to the time derivative of the diagonal 
density matrix elements in the following way: 

(^ L )dias =r»(nf -pu), (19) 

two terms corresponding to generation and absorption of the plasmons. The temperature of Bose distribution function 
nf corresponds to the electron temperature of the ith layer. Non-diagonal matrix elements aquire a damping equally 
from both layers 



The system of equations that incorporates both dissipation and resonant tunneling reads as follows: 



Pij (20) 
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9pn r , B \,- SuJ r \ 
= r^ri! - pii) + i—(pi2 - Pn) 

= T 2 {n 2 - p 22 ) + «^-(P2i - P12) 

dpu (ri + r 2 ) ,.&*). , 

~m= 2 — ^ + iT (pn-p 22 ) 

where P21 — p*2- The stationary solution takes the form 

b T 2 6w 2 (nf - ng) 
Pu = % - — (21) 

1 1 + I 2 + 1 il 2 

Expression for P22 is obtained by reverting indices 1,2. 

We can now evaluate the drag force acting on electrons of each layer by equating it to a momentum flow between 
the layers. We sum over modes with all possible q and obtain 

f = - V nA diss = - y ft q FlF2 ^1 (22) 

where T, n B , 5ui may be q-dependent. 

We assume that the current flows in the layer 2 and the drag force in the layer 1 is equilibrated by the electric field. 
The transrestivity is essentially the ratio of this field to that current. The effect of the current is that the n B is the 
equilibrium Bose distribution in the reference frame where the electrons of the second layer are in average at rest, 
rather than in the laboratory reference frame. So that 

ni(e{q)) = /s(e(q) - h(v d „ ft q)) / B (e(q)) - h(v drift q)-^- (23) 

Vdrift being the drift velocity. 

Substituting ( p3|) in (^2|) we obtain 

- h v dnf i] Ti+T2Slu 2 +TiT2 d£ ( 24 ) 

The last things to note are that F = en\E and I = eri2Vdrift- If we use this and set ri = T 2 — T we reproduce 
Eq.©. 



VI. RESULTS: HIGH TEMPERATURES 



In this section we consider the drag resistance at temperatures T ^> Tllu c , which are sufficiently high to excite 
magnetoplasmons and electrons in many Landau levels. As we can see in Fig. 1, at high temperature we encounter 
at least three distinct regions with different temperature and magnetic field exponents. 

First of all, we shall explain why there are so many regions. If we compare imaginary part of polarization function 
with and without magnetic field, we see that their values averaged over frequency intervals bigger than Tiuj c are the 
same < (Im%) >=< (Imx(H = 0)) >. However, Im% = beyond narrow adsorption bands (Fig. 3). This means that 
whinin the bands Im% is significantly enhanced in comparison with its zero-field value, typically by a factor u c /A. 
Surprisingly enough, the enhancement of Im% can lead both to enhancement and suppresion of the drag. 

If Imx <C v, Rex ~ v. The denominator £ in Eq. [I] which is responsible for screening of interlayer potential 
is the same as without magnetic field. We refer to this situation as to normal screening regime. In this regime, 
the transresistance is enhanced in comparison to its value without magnetic field, since the effect is proportional to 
< (Imx) 2 »< (Imx) > 2 ~< (Imx(-ff = 0)) > 2 . 

At further increase of Imx, R- e X develops as well so that both | Rex |, Imx become bigger than v. The denominator 
£ strongly increases. That efficiently screens out the inter-layer interaction and leads to drastic descrease of the 
transresistance. That we will call over screening. 

However, £ can also decrease with increasing Imx an( I pass zero. Near this line, the inter-layer interaction is greatly 
increased. This is where the magnetoplasmon contribution dominates. 

The actual value of the drag effect is thus determined by interplay of these three competing tendencies. 
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Let us first evaluate the magnetoplasmon contribution. Substituting expressions ([i"3|), (14) to Eq.(|l7|) we notice 
that the integrand has a sharp extremum near n sa so we can formally integrate over k in infinite limits. This yields 
the relation which is valid for all regions, 



w» = -0 0Q221 hA ^ a 4 V {h ^ /T)4 ( a ^ ,T X 3 " 

Pl2 e 2 n in2 d*T\d) ^su& 2 {hw j /2T)\sx a h{ahu j /T) j 



(25) 



where a ~ 0.19(Tw c /A£;o)- 

This expression can be simplified further. The region I (huj c « T « EoA/huj c ) corresponds to a <C 1. Since 



7ia; c <C T the characteristic values of j in the sum (25) are much bigger than unity. Therefore we can convert the sum 
over j into the integral. Since f °° dxx 4 sinh _2 (iz;/2) = 167r 4 /15, we obtain 



-0.0003^^ (^) 3 ' S . (26, 
e i n\n 2 a i A^Eq V a / 



The region II is defined by inequalities T ^> EqA/ui c ^> w c . Here a ^> 1. As a result, the characteristic frequencies 
here Qj — EqA/lo c <C T (see Eq.(^5])). These frequencies, however, are still much larger than cyclotron frequency: 
u)j S> w c . This enables us again to introduce the continuous variable x = othujj/T and convert the sum into the 
integral: 



In region III (u> c ^> y/EoS) the j — 1 term dominates the sum. The corresponding expression for the plasmon drag 
resistance is exponentially small, 

pg * -2.7 x 10- 6 2 ^ ( ^ e) ^ 2 (^) 3/2 exp(-0.28(^ c ) 2 /Ai?o). (28) 
e' L n\n 2 dr A 9 / 2 Eq \ a J 

This is due to the fact that the value of q which is needed to bring the magnetoplasmon pole to the vicinity of the 
Landau level is large compared to 1/d. Similar situation occurs in the region IV, with the exponential suppression 
being due to low temperature, 

* * - 112 x 10 ~ 5 ^rb§f (^r^(-^m- <*> 

We will see below that even the exponentially suppressed magnetoplasmon contribution can efficiently compete 
with the quasiparticle one. 

Now we will estimate quasiparticle contribution in all three regions. 

Let us first consider region I (uj c «T« EqA/lo c ). In this parameter interval the characteristic values of Tvjj ~ T 
and qd ~ 1. It follows from Eqs.0,|| that Rex = v and Im% <C v. Thus we are in the normal screening regime. Using 
the fact that V(q)v 3> 1 (in other terms, qas <C 1) we obtain from Eq.|]: 

^--0.007254 1 M IT- (30) 

This magnetotransresistance is bigger than the zero-field value (Eq.||) by a factor of lu c /A 1. This is due to 
the discreteness of the electron spectrum in the magnetic field when the density of states within the Landau level 
remarkably increases. 

Region II (ui c -C E^A/lo c <C T). Though the main contribution to the drag in this region is due to the magnetoplas- 
mon mechanism, it is instructive to give an estimation of the quasiparticle contribution. As in the region I, q ~ 1/d. 
The characteristic frequency, however, is restricted by the value U) C ut ~ EqA/lo c which is much smaller than the tem- 
perature (though still much larger than the cyclotron frequency). The reason is the overscreening. For estimations, 
we approximate Imx near the cyclotron resonance by its value from Eq.0 at qd ~ 1. That gives Imx ~ vui/uicut- 
A good estimation for Rex near the resonace is v for v 3> Imx and Imx otherwise. Thus the integrand in Eq. ([!]), 
(Imx) 2 / | £ | 2 cx (Imx) 2 / (Re 2 x + Im 2 x) 2 achieves a maximum at uj cu ±. Since the interval of integration over fre- 
quency is effectively smaller than the temperature, the transresistance exhibits a linear temperature dependence: 
pll ~ (?i/e 2 )(l/n 1 n 2 d 4 )(a 2 J /d 2 )(r/i? ). 
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Let us compare this with B = case. There, the linear temperature dependence startsa at T > Eq. This is because 
the absorbtion of energy due to the Landau damping mechanism is possible only for frequencies smaller than vf/cI. 
Thermal frequencies are ineffective, because corresponding phase velocities are larger than vf- In the presence of a 
magnetic field, the role of a cut-off energy is taken by a much smaller value EqA/u c . So that, the linear temperature 
dependence of the transresistance starts at much lower temperatures. 

Region III (\/EqA <C fvu c <C T). Here u cut becomes smaller than u c . There is overscreening at qd ~ 
1 for all cyclotron resonances. The contribution of resonances to integral Eq.(Q) is of the order of pff ~ 
(^/e 2 )(l/n 1 n 2 d 4 )(i? c 2 a 2 B /d 4 )(TA 3 / Wc 4 ). 

The main contribution is determined by the frequencies u ~ A <C u c . So that the quasiparticles are created within 
the same Landau level. For these frequencies we always have the normal screening situation. We can set Rex = v. 
The Imx is determined by the j = term of Eq.(j?]). We obtain 

= _h^J_ f (Pqq 2 1 Ta% f 2A du , us ^ K 1 u\Ta\ 

Pl2 e 2 ir* ni n 2 J (2tt) 2 s inh 2 (qd) A?R cl R c2 J 2ir m \2AJ- ' e 2 mn 2 d 4 Av F1 v F 2 1 ' 

It is interesting to note that this magnetic and temperature dependence coincides precisely with the observed one@ if 
we assume that A does not depend on magnetic field. Indeed, the magnetic dependence of A is rather weak. These 
experiments were performed in rather strong magnetic field where only few Landau levels were occupied so that one 
should not expect quantitative agreement with our calculations. From the other hand, the linear T-dependence is 
rather remarkable. We believe that in any case it indicates a reduction of typical excitation energies to values much 
smaller than T, possibly due to overscreening at energies of the order of T. 

Now we are in position to compare quasiparticle and magnetoplasmon contribution and thus to set the bor- 
ders of the gray-shaded regions in Fig. 1. The magnetoplasmon contribution dominates throughout the region 
II. In the region I we compare expressions ( p6| ) and (|30|). Magnetoplasmon contribution dominates provided 
T ~ EoA/hu c (aB j 'd) 1 / 4 . Since experimentally d ~ as this happens in fact close to the border between the re- 
gions I and II T ~ EoA/huj c . The exponentially small magnetoplasmon contribution given by ( p8| ) competes with 
( pl| ) in the region III and dominates provided 0.28w 2 /(Ai?o) < ln((w 2 /i?oA) 7 / 2 (d/aB) 1/ ' 2 ). It also dominates in the 

region IV if T/huj c > 1/ ln(w 4 Td 1 / 2 /(i? 2 A 3 a B 2 )). The latter condition is obtained by comparing expressions ( |29| ) 
and @. 

VII. RESULTS FOR LOW TEMPERATURES 

In this section we present our results for low temperatures T -C Tiu c . Owing to energy limitations, only the states of 
the upper partially filled Landau level arc involved into the drag. This makes the transresistance sensitive to concrete 
value of this filling factor. In addition to slow dependence on magnetic field, the drag effect exhibits oscillatory 
dependence on inverse magnetic field related to the filling factor. The detailed study of these Shubnikov-de Haas 
oscillations is beyond the limits of the present work and will be presented elsewhere. This is why we provide here 
analythical results for the regions IV and VI only. As to the regions V and VII, we present below estimations of the 
transresistance for typical filling factors rather than detailed analythical results. 

In the region IV (A <C T fiu c <C \JTEq) the situation is the most straightforward one since it follows from 
Eq-@ that Imx ^ v - Thus we encounter here the normal screening and we may set Rex = v. We obtain for the 
transresistance: 

PYf = -0.0011 A 1 a e - /„J/„ 2 (1 - /„ 2 ) (32) 

f ni 2 being filling factors in the layers. 



Even if T 3> A, we encounter overscreening in the region V (\/TEq <C tiu c <C Eq). For estimations, we set 
Imx ~ R- e X ^ v an d obtain 

P12 - (?»/e 2 )(l/7i 1 7i 2 d 4 )(a 2 3 i? cl i? c2 /d 4 )(r 3 /^ c 2 A). (33) 

As a consequence of the overscreening, the transresistance rapidly decreases with increasing magnetic field. 

At low temperature T <C A the integral in Eq. 1 is contributed by u ~ T. This gives rise to featureless T 2 
temperature dependence of Eq. || with the coefficient depending on magnetic field. The region VI (T <C A, lo c <C 
\/EqA) again corresponds to normal screening. We take Im from Eq.dTtl) and set Rex = v. This yields 
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d xx - 0031 ^ 1 °* r2 ^ fl (Ml ~ C " )2 Vl ( ^ 2 ~ £ " ) % (U) 

Pl2 - U -°^ 1 e 2 nin2d 4 WF1 ^ 2 A 4 C 1 A 2 )(1 A 2 ) ( 34 ) 

where /ii.2 are chemical potentials in the layers. They are related to filling factors f ni ,n 2 by means of 



/„ = 1/2 + (!/*)[(! - (Ml A2 £ " )2 )yi-(l- (Afl A2 £ " )2 ) 2 + arcsin(l - ^ J n? )]. (35) 

For a typical filling factor, the effect is bigger than zero filed transresistance by a factor of (w c /A) 4 . 

With increasing the magnetic field, we enter the region of overscreening (VII, T <C A, uj c ^> ^EqA.) Again we set 
Rex ^> v and obtain 

Pi2 ~ (Ve 2 )(l/nin2d 4 )(4i? c ii?c 2 /d 4 )(T 2 K 2 ). (36) 
that decreases with increasing magnetic field. 
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APPENDIX: 

In this work we disregard rapidly oscillating parts of polarization function given by Eq. |. This means we approx- 
imate Bessel function squares in Eqj^ in the limit of qR c 3> 1 as J 2 n (qR c ) « l/nqR c rather than as 

JU<l R c) « -4t(1 + cos(2gi? c - Trm - J)), (Al) 

which is the mathematically correct expression. Let us explain why. 

First let us note that if we take these oscillating parts into account it would significantly alter our results. In the 
normal screening regime it would give an extra factor of 3/2, since the answer is proportional to (Im%) 2 . The answer 
would change even more drastically in the overscreening regime. The point is that oscillating terms would set the 
polarization function to (almost) zero for q corresponding to zeros of the Bessel function. No overscreening would 
occur near these points and their close vicinity would dominate the drag. 

All this would lead to very sophisticated and extremely unstable picture of the drag effect. Fortunately, we are able 
to present some arguments that allow to disregard oscillating terms in polarization function. 

The physical origin of the oscillating terms can be best understood in the language of semiclassical electron trajec- 
tories in magnetic field. A classical trajectrory can not move from the starting point further than 2R C . This is why 
the polarization function in coordinate representation, x( x i x> ) has a sharp edge at | x — x' |= 2R C . This gives rise to 
Fourier components ~ cos(2qR c ). 

If the edge is not sharp, the oscillating part is exponentially suppressed. The suppression is of the order of 
exp(— q 2 (SR) 2 ), 5R being a typical rounding of the edge. By virtue SR is the typical uncertainty of the coordinate 
of the electron that makes a half of Larmor circle. Such an uncertainty can be of quantum-mechanical origin. In 
this case we estimate SR ~ \jR c /kp *C R c - Another cause of uncertaincy may be small- angle scattering by smooth 
potential fluctuations in the heterostructure. For this case we estimate SR ~ y/R^/l sa , where transport mean free 
path l sa R c . It is interesting to note that scattering on point- like defects does not contribute to SR provided 

LU C T » 1. 

Now we note that typical q contributing to the drag resistance are of the order of 1/d. We conclude that the 
oscillating part is exponentially suppressed provided d < y/ R c /kF or d < R C ^J R c /lsa- We assume that at least one 
of these conditions is fulfilled. 
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FIG. 1. Seven regions of different analytical behaviour of the transresistance in the intermediate magnetic field regime. 
Note log- log scale. First number in each region corresponds to temperature exponent, second number indicates magnetic field 
exponent. Vertical line corresponds to the condition (huj c ) 2 ~ EqA. 



FIG. 2. The functions Xim(x) and X r e{x) determine the shape of Imx, Rex in the vicinity of each cyclotrone resonance. 

FIG. 3. Electron adsorption in u),q plane. Left plane: electron adsorption occurs i. in narrow strips around to = and 
cyclotron resonances (particle-hole continium) ii. on the magnetoplasmon dispersion curves. Right plane: intersection of the 
magnetoplasmon dispersion curve and the edge of the strip at a smaller scale. We illustrate splitting of the magnetoplasmons 
and the level width Y they acquire inside the strip. 
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